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Abstract 

We study the possibility to analyse the data on inclusive prompt photon production at 
first LHC runs in the framework of the fc^-factorization QCD approach. Our consideration is 
based on the amplitude for the production of a single photon associated with a quark pair in 
the fusion of two off-shell gluons. The quark component is taken into account separately using 
the quark-gluon Compton scatterring and quark-antiquark annihilation QCD subprocesses. 
The unintegrated parton densities in a proton are determined using the Kimber-Martin- 
Ryskin (KMR) prescription as well as the CCFM evolution equation. 

PACS number(s): 12.38.-t, 13.85.-t 
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The present note is motivated by the very recent measurements [IJ[2] performed by the 
CMS and ATLAS collaborations at the LHC where the isolated photon cross sections have 
been presented as functions of their transverse momentum at y/s = 7 TeV for the first time. 

The production of prompt photon^] in hadron-hadron collisions at high energies is a 
subject of pointed discussions up to now. The theoretical and experimental investigations of 
such processes provide a direct probe of the hard subprocess dynamics since the produced 
photons are largely insensitive to the effects of final-state hadronization. The cross sections 
of these processes are strongly sensitive to the parton (quark and gluon) content of a proton 
since isolated prompt photons can be produced mainly via quark-gluon Compton scattering 
or quark- ant iquark annihilation (at LO). 

We note, however, that a complete theoretical description of the Tevatron data within the 
QCD is an open question still (see j3jll] and references therein). It was demonstrated [5j|6] 
that the overall description of these data in the framework of the standard QCD can be 
achieved by introducing some additional intrinsic transverse momentum of incoming 
partons, which is usually assumed to have a Gaussian-like distribution. The average value 
of this k? increases from kx ~ 1 GeV to more than kx ~ 3 GeV in hard scattering processes 
as the y/s increases from UA6 to Tevatron energies [51E]- The importance of including the 
gluon emission through the resummation formalism was recognized and only recently this 
approach has been developed for inclusive prompt photon production [8HTT]. 

In the framework of the /^-factorization approach of QCD [12], which is of primary con- 
sideration in this note, the transverse momentum of incoming partons occurs in a natural 
way. In this approach, the non-zero partonic kx is generated perturbatively in the course 
of non-collinear parton evolution via the corresponding (usually BFKL [13J or CCFM |14j ) 
evolution equations. A detailed description of the /c^-factorization can be found, for exam- 
ple, in reviews [15J. Studies of the prompt photon production at hadronic colliders in the 
framework of hr- factorization approach have been done in many papers [T6H20] . So, investi- 
gations [T6HT8] have been based on the leading-order matrix elements of quark-gluon Comp- 
ton scattering and quark- ant iquark annihilation subprocesses. An important component of 
these calculations is the unintegrated {kx- dependent) quark distributions in a proton. At 
present, these densities are available in the framework of Kimber-Martin-Ryskin (KMR) [21] 
approach only since there are some theoretical difficulties in obtaining the quark distribu- 
tions immediately from CCFM or BFKL equations (see, for example, [15] for more details). 
In our previous investigations [EJEn] we tried an another way. The main idea of [T91I20] 
was in reexpressing the quark contributions in terms of gluon ones using the higher-order 
off-shell gluon-gluon fusion matrix elements, namely g* + g* — > 7 + q + q. Thus, we reduced 
the problem of poorly known and poorly calculable unintegrated quark distributions to much 
better investigated gluon densities. The corresponding contributions from the valence quarks 
have been taken into account separately. We obtained [IHJEO] a reasonably good agreement 
between the /^-factorization predictions and the Tevatron data on the inclusive prompt 
photon production cross cross sections in both central and forward photon pseudo-rapidity 
regions. Based on the results [T9|l20]. in the present note we give the first systematic analysis 
of data [TJ[2] taken by the CMS and ATLAS collaborations at the LHC in the framework 
of the /^-factorization approach. Moreover, we improve our previous predictions [121120] by 

1 Usually the photons are called prompt if they are coupled to the interacting quarks. 
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taking into account the transverse momentum of incoming off-shell quarks in a proper way. 
In other aspects we will strictly follow the approach described in jT^I2U]. 

We only briefly recall here the corner-stones of proposed theoretical scheme. The starting 
point of consideration is the usual leading order 0(a) and 0(aa s ) off-shell subprocesses, 
namely q* + g* — > 7 + q and q* + q* — > 7 + g. These subprocesses are strongly depend on the 
unintegrated quark densities in a proton, f q ( ). The unintegrated quark distributions 

include the ones of valence quarks ff^( ), sea quarks appearing at the last step of 

the gluon evolution f^(x, k 2 -,, /i 2 ) and sea quarks coming from the earlier (second-to-last, 
third-to-last and other) gluon splittings f^ s '( X j k<2~> ; jJi ). In the proposed approach [191 120] 
we simulate the last gluon splitting by the higher-order 0(aa 2 ) off-shell matrix elements 
g*+g* — > 7+g+g. In this way we take into account the contributions from the f^ 9 \ 
To estimate the contributions from the ff} v \x, k^, fj?) and ff\x, k 2 -., fi 2 ) we use the specific 
properties of the KMR formalism which enables us to discriminate between the various 
components of the unintegrated quark densities (see below). 

Thus, the proposed scheme results to the following partonic subprocesses^]: 

9* + 9* + q + (1) 

q' M + 9*^l + q, (2) 
q <v,s) + y*(v,s) ^ 1 + g , (3) 

To be precise, the gluon-gluon fusion (1) takes into account the contribution of the q*^ + 
q*^ annihilation subprocesse, and the valence and sea quark-gluon scattering (2) take into 
account the q*( v > + q*^ and g*W 4- q*W mechanisms. Evaluation of the off-shell matrix 
elements of subprocesses (1) — (3) is straightforward and the analytical expressions have 
been listed in [T9]. Here we only would like to mention two technical points. First, according 
to the fc^-factorization prescription [12], the summation over the incoming off-shell gluon 
polarizations in (1) and (2) is carried with ^e^e v = k^k^/k 2 -., where kr is the gluon 
transverse momentum. Second, when we calculate the matrix element squared, the spin 
density matrix for off-shell spinors in (2) and (3) is taken in the form u(p)u(p) = xp p [22J, 
where x is the fraction of initial proton longitudinal momentum p p . In all other respects our 
calculations follow the standard Feynman rules. Since the expression for the off-shell quark 
spin density matrix has been derived in the massless approximation, numerically we neglect 
the charmed quark mass. 

According to the /^-factorization theorem, to calculate the cross section of the prompt 
photon production one should convolute the off-shell partonic cross sections (1) — (3) with 
the relevant unintegrated quark and/or gluon distributions in a proton: 

o = J °ij( x ii x 2, k 2 T , k 2 2T ) fi(xi, k 2 T , n 2 )fj(x 2 , \4t, Z^) dxidx 2 dk 2 1T dkl T , (4) 

i,3=q,9 

where a*j(xi, x 2 , k 2 T , k| T ) is the relevant partonic cross section. The initial off-shell partons 
have fractions X\ and x 2 of initial protons longitudinal momenta and non-zero transverse 



2 We will neglect the contributions from the so-called fragmentation mechanisms. It is because after 
applying the isolation cut (see [i|2] ) these contributions amount only to about 10% of the visible cross section. 
The isolation requirement and additional conditions which preserve our calculations from divergences have 
been specially discussed in [TW20] . 
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momenta and k 2 ^. The analytical expressions for the contributions of subprocesses (1) 
— (3) are given in [T9] . 

Concerning the unintegrated parton densities, we use the CCFM-evolved gluon and va- 
lence quark distributions derived in [23] and [24] , respectively. To determine unintegrated 
parton densities we also apply the KMR approximation [21]. The KMR approach is the 
formalism to construct the unintegrated parton distributions f a ( ) from the known 

conventional parton distributions xa(x,fi 2 ), where a = g or a = q. In this approximation, 
the unintegrated quark and gluon distributions are given by [2T] 




where P a b(z) are the usual unregulated LO DGLAP splitting functions. The theta functions 
which appear in (5) and (6) imply the angular-ordering constraint A = + (krl) specif- 
ically to the last evolution step to regulate the soft gluon singularities. For other evolution 
steps, the strong ordering in transverse momentum within the DGLAP equations automat- 
ically ensures angular ordering^. The Sudakov form factors T g (k^,/i 2 ) and T 9 (k^,/x 2 ) which 
appears in (5) and (6) enable us to include logarithmic loop corrections to the calculated 
cross sections. 

Note that the function f q (x, k 2 ^, ll 2 ) in (5) represents the total quark distribution function 
in a proton. Modifying (5) in such a way that only the first term is kept and the second term 
omitted, we switch the last gluon splitting off, thus excluding the fjf\x, k 2 -,, ll 2 ) component. 
Taking the difference between the quark and antiquark densities we extract the valence quark 
component fj^ v \x, ky, /x 2 ) = f q (x, k 2 -, li 2 ) — fg(x, k 2 -, fx 2 ). Finally, keeping only sea quark in 
first term of (5) we remove the valence quarks from the evolution ladder. In this way only 
the f( s \x, ky, li 2 ) contributions to the f q (x, k^, li 2 ) are taken into account. 

Other essential parameters were taken as follows: renormalization and factorization scales 
ll = ^Ej< (where we vary the parameter £ between 1/2 and 2 about the default value £ = 1 
in order to estimate the scale uncertainties of our calculations), LO formula for the strong 
coupling constant a s (/i 2 ) with rif = 4 massless quark flavours and Aq CD = 200 MeV, such 
that a s (M|) = 0.1232. The multidimensional integration in (4) has been performed by 
means of the Monte Carlo technique, using the routine VEGAS [26]. The full C++ code 
is available from the authors on request^. This code is practically identical to that used 
in [19l[20] with exception that now we apply it to the LHC conditions. 

The results of our calculations are shown in Fig. 1 where we confront the calculated 
differential cross section of the inclusive isolated prompt photon production (as a function 

3 Numerically, in (5) and (6) we have applied the Gliick-Reya-Vogt (GRV) parton distributions [25] . 
4 lipatov@ theory, sinp . msu . ru 
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of the photon transverse energy E T ) with the first LHC data [HE] taken by the CMS and 
ATLAS collaborations at y/s = 7 TeV. In Fig. 2 we show the data/theory ratio of our 
predictions. The CMS data refer to the kinematic region defined by E T > 21 GeV and 
\rp\ < 1.45. The ATLAS data have been measured at 15 < E} < 100 GeV, < 0.6, 
0.6 < |y 7 | < 1.37 and 1.52 < |y 7 | < 1.81, respectively. In Figs. 1 and 2, the solid and 
dash-dotted histograms are obtained with the KMR and CCFM parton densities by fixing 
both the factorization and normalization scales at the default value \x = E T . The upper 
and lower dashed histograms correspond to the scale variation in KMR predictions as it was 
described above. One can see that the predictions obtained with the KMR partons tend to 
slightly overestimate the measured cross sections, but they are coincide with the data within 
the theoretical uncertainties. Contrary, predictions based on the CCFM evolution lie below 
the LHC data in a wide E T range. The main reason of such behaviour is connected with 
taking into account the contributions from the sea quarks originating from the earlier steps 
of the evolution cascade in the former case. However, to avoid the possible double counting 
these contributions are not taken into account in the CCFM calculations since part of them 
can be already included into the CCFM results (via initial parton distributions which enter 
to the CCFM equation). Note that the similar effect has been observed [22] in the prompt 
photon photoproduction at HERA. 

The relative contributions of subprocesses (1) — (3) to the prompt photon cross section 
are shown in Fig. 3. Here we use the KMR parton densities for illustration. One can see that 
at low Et the gluon fusion subprocess gives the dominant contribution whereas the QCD 
Compton and quark-antiquark annihilation subprocesses are mostly significant at higher Et- 
The contribution from the quarks appearing from the earlier evolution steps, f^ s '(x, k T , /i 2 ), 
is also significant in both central and more forward (as it was measured by the ATLAS 
collaboration) rapidities. 

To conclude, in the present note we apply the CCFM and KMR unintegrated parton den- 
sities to the analysis of the first experimental data on the prompt photon production in pp 
collisions taken by the CMS and ATLAS collaborations at the LHC. We have obtained a rea- 
sonably good agreement between our predictions and the LHC data and have demonstrated 
that the quarks are important even at the LHC energies and therefore should be properly 
included into the non-collinear evolution equations. It is important for further studies of 
small-x physics at the LHC, in particular, for investigation of Drell-Yan pair production 
which give us possibilities to test the region of very low x (up to x ~ 10~ 5 ). 
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Figure 1: The differential cross sections of the inclusive prompt photon production in pp 
collisions as a function of calculated at y/s = 7 TeV. The solid histogram corresponds 
to the KMR predictions at the default scale /i = Ej,, whereas the upper and lower dashed 
histograms correspond to scale variations described in the text. The dash-dotted histogram 
corresponds to the predictions obtained with the CCFM parton densities. The experimental 
data are from CMS [1] and ATLAS 0. 
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Figure 2: The theory/data ratio of the inclusive prompt photon production in pp collisions 
as a function of Ej, calculated at y/s = 7 TeV. Notation of the histograms is the same as in 
Fig. 1. The shaded band represent the experimental uncertainties. 
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Figure 3: Different contributions to the cross section of the inclusive prompt photon produc- 
tion in pp collisions as a function of Ej, calculated at y/s = 7 TeV. The dashed, short dashed 
and dotted histograms correspond to the contributions from the q*g* — > 97, g*g* — > qq^ 
and q*q* — > subprocesses, respectively. The dash dotted histogram corresponds to the 
fq S \x, k^, /i 2 ) contribution. The solid histogram represents the sum of all these components. 
We use the KMR parton densities for illustration. The experimental data are from CMS [1] 
and ATLAS [2]. 
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